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Abstract 

We propose the gravitino dark matter in the gravity mediated supersymmetry breaking sce¬ 
nario. The mass hierarchies between the gravitino and other superparticles can be achieved 
by the nontrivial Kahler metric of the SUSY breaking field. As a concrete model, we con¬ 
sider the five-dimensional supergravity model in which moduli are stabilized, and then one 
of the moduli induces the slow-roll inflation. It is found that the relic abundance of gravitino 
and the Higgs boson mass reside in the allowed range without a severe fine-tuning. 
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1 Introduction 


The low-scale supersymmetry (SUSY) is an attractive scenario which not only protects the 
mass of the Higgs boson from the large radiative corrections but also gives the dark matter 
candidates. In addition to it, the existence of supersymmetry is also motivated in the string 
theory which is expected as the ultraviolet completions of the standard model (SM). This is 
because the SUSY guarantees the absence of tachyons in the string theory. 

In the minimal supersymmetric standard model (MSSM), the large radiative corrections are 
indicated by the observed Higgs boson mass within ranges between 124.4 and 126.8 GeV [Tj. One 
of the solutions to raise the Higgs boson mass in the MSSM is the high-scale SUSY-breaking 
scenario, and then the SUSY flavor and CP problems can be also solved at the same time. 
However, this scenario brings the tuning problem to the MSSM in order to realize the successful 
electroweak (EW) symmetry breaking. By contrast, there is another solution to raise the Higgs 
boson mass by the nature of left-right mixing of the top squarks. As pointed out in Ref. [2], the 
nonuniversal gaugino masses at the grand unification theory (GUT) scale Mgut — 2 x 10 16 GeV 
lead to the maximal mixing of the top squarks and then, the realistic Higgs boson mass can be 
achieved without a severe fine-tuning by the structure of the renormalization group equations 
in the MSSM. Throughout this paper, we focus on this low-scale SUSY breaking scenario. 

The SUSY-breaking scenarios are mainly categorized into the gravity mediation 0 , gauge 
mediation [4J and the anomaly mediation [5], For any mediation mechanisms, the gravitino 
mass is sensitive to the cosmological problem, e.g., the cosmological gravitino problem [6]. If 
the gravitino is not stable, the mass of the gravitino should be larger than 0(10 TeV) in order 
to be consistent with the successful big bang nucleosynthesis (BBN). The lower limit of the 
gravitino mass depends on the reheating temperature; for more details see Refs. PIS El El Eo] Q 
Therefore, before discussing our considered situation, we comment on several SUSY-breaking 
scenarios, focusing on the mass of the gravitino. 

In the gauge mediated SUSY-breaking scenario, the dark matter candidate is the ultralight 
gravitino of mass 7713/2 *C 0(1 GeV) under the low-scale SUSY-breaking. Note that if the 
gravitino mass is larger than this scale, it is expected that the gravitational interactions give 
the sizable effects to the dynamics of the SUSY-breaking sector as well as the visible sector. 
In the pure anomaly mediated SUSY-breaking scenario, the wino-like neutralino is likely to 
be the dark matter candidate due to the structure of the beta functions in the MSSM [12]. 
However, the recent results of the LHC experiments [13j indicate the TeV scale gluino mass; in 
other words, the large mass of the gravitino m 3 / 2 — 0(100 TeV) is required in the framework 
of anomaly mediation. In the mirage mediation [14] . the mixed neutralino would be the dark 
matter candidate and the large gravitino mass above 0(10 TeV) is expected in the light of the 
cosmological gravitino problem. In the gravity mediation, the neutralino dark matter is often 
considered under the large gravitino mass above 0(10 TeV) with high-scale SUSY breaking, 
otherwise, SUSY flavor violations arise due to the flavor dependent interactions. 

In this paper, we consider the gravity mediated SUSY-breaking scenario which is compatible 
with the low-scale SUSY and observed Higgs boson mass without the cosmological gravitino 
and SUSY flavor problems. In general, the gravity mediation connects the scale of the grav- 

1 It is also possible to consider the light gravitino in the extension of the MSSM. See, e.g., Ref. m- 
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itino mass with that of the supersymmetric particle, because the origin of soft SUSY-breaking 


terms is only the gravitational interactions. Therefore, it seems to be difficult to solve the 
cosmological gravitino problem with the low-scale SUSY-breaking scenario. In order to realize 
the low-scale SUSY without the cosmological gravitino problem, we propose the mechanism 
to generate the mass hierarchies between the gravitino and the other sparticles based on the 
framework of a four-dimensional Af — 1 supergravity (4D Af = 1 SUGRA). Especially, we focus 
on the case that the gravitino is the lightest supersymmetric particle (LSP) whose mass is of 
0(100 GeV). Since such a stable gravitino is much heavier than that predicted by the gauge 
mediated SUSY-breaking scenario, this would be the typical feature of the gravity mediation. 
There are some studies for the gravitino dark matter with assumed sparticle spectra that focus 
on the cosmological implications and it is then found that the next-to-the-lightest supersym¬ 
metric particle (NLSP) is severely constrained. (See, e.g., Refs. [TUJ [15, [16] H7].) In order to 
determine the relevant higher-dimensional operators in 4D Af = 1 SUGRA, we consider a five¬ 
dimensional supergravity (5D SUGRA) compactified on an orbifold S 1 /Z 2 . In the framework 
of 5D SUGRA, the successful inflation mechanism as well as the moduli stabilization can be 
realized as suggested in Refs. [18, [19]. The dynamics of inflaton and moduli are important to 
evaluate the abundance of the gravitino produced via the inflaton and moduli decay into the 
gravitino. Furthermore, the Yukawa hierarchies of elementary particles can be realized without 
a severe fine-tuning by employing the localized wavefunction of quarks, leptons and Higgs in 
the fifth dimension ra. 

The following sections are organized as follows. In Sec. El we discuss how to realize the mass 
hierarchies between the gravitino and other supersymmetric particles in 4D Af = 1 SUGRA. As 
a concrete model, in Sec. [3] we briefly review the structure of 5D SUGRA on S 1 /Z 2 and then 
the gravitino can be the dark matter candidate. Thanks to the detailed moduli stabilization as 
well as the inflation mechanism, one can discuss the nonthermal productions of the gravitino 
via the moduli and inflaton decay after the inflation. After that, in Sec. 0 we evaluate the relic 
abundance of the gravitino and the Higgs boson mass with a severe fine-tuning. The obtained 
results are consistent with the cosmological observations as well as the collider experiments. 
Finally, Appendices [A] and [B] denote the detailed derivations of the scalar potential around the 
vacuum and during the inflation, respectively. 

2 The mass hierarchies between the gravitino and other 
sparticles 

In this section, we show how to realize the mass hierarchies between the gravitino and other 
sparticles in the framework of 4D Af = 1 SUGRA. The scalar potential in 4D Af = 1 SUGRA 
is given by 



K I jF I F J -?,e K/2 \W\\ 


( 1 ) 
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where K and W are the Kahler and superpotential, respectively. DjW = Wj + KjW with 
Wj = dW/dQ 1 , Kj = dK/dQ 1 are the Kahler covariant derivatives of the superpotential for 
the scalar components of the chiral superfields Q 1 , F 1 = —e K ^ 2 K IJ DjW are the F-terms of Q 1 
and K IJ are the inverse of Kahler metric Kjj = djdjK. Here and hereafter, we set the Planck 
unit Mpi = 1, unless we specify it. The vanishing cosmological constant (V) = 0 is rewritten 
in the following form: 


"> 3/2 = \{KxxF x F x ), (2) 

where 7713/2 = e < ' K ' > ^ 2 (W) is the gravitino mass. It is then assumed that the SUSY is broken 
by the single chiral superfield whereas the soft SUSY-breaking masses of the gauginos and 
scalar components of the chiral superfields Q 1 are given by 

M a = (F x d x \n(Ref a )), 

ml = -(F x F x d x d x lnY <p r), (3) 

where f a , a = U(l)y, SU(2) L , SU(3)c are the gauge kinetic functions of the standard model 
gauge groups whose vacuum expectation values (VEVs) determine the size of gauge couplings. 
Yqi are some nontrivial functions for the kinetic term of Q 1 which can be severely constrained 
by the flavor structure of elementary particles as can be seen later. From the above Eqs. (|2j) 
and (f3l) . the nontrivial Kahler metric of the SUSY-breaking field X gives rise to two nontrivial 
possibilities: 

• The gravitino dark matter: 

I 11 the case of (K xx ) 1, the gravitino mass is smaller than the soft SUSY-breaking 
masses for any value of the F-term (F x ). Then it is possible to consider the gravitino 
dark matter in the gravity mediated SUSY-breaking scenario with TeV scale gauginos and 
sparticles. It is then assumed that the derivatives of the gauge kinetic function d x Kef a 
and the kinetic term of Q 1 , d x d x lnY^i satisfy the certain conditions in order to obtain 
the gravitino dark matter. This is because the renormalization group (RG) effects are 
significant to discuss the sparticle spectrum. Such conditions are discussed in the case of 
constrained MSSM (CMSSM) [2Tj . 

The stable gravitino would be consistent with the thermal history of the universe, even 
if the decays of NLSP do not spoil the success of BBN [TDJ [151 'US EJ and at the same 
time, the relic abundance of the gravitino should not be larger than that reported by 
the Planck Collaboration [22]. In any case, the stable gravitino is favored in the light 
of naturalness, because the F-term of the SUSY-breaking field can be taken as a usual 
low-scale SUSY-breaking scenario which soften the divergences for the Higgs boson mass. 
Note that the small Kahler metric of the held X should be ensured in order not to be 
below that generated by the loop and/or higher derivative corrections. 

2 It is straightforward to extend our situation in multiple SUSY-breaking fields. 
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• Other dark matter candidates: 

By contrast, in the case of (K xx ) G if i s expected that the gravitino is heavier than the 
other sparticles for any value of the F-term ( F x ) with d x Ref a and d x d x \nY$i of order 
unity. Thus, one can solve the cosmological gravitino problem with the low-scale SUSY- 
breaking scenario. Then, the gravitino mass can be chosen as above 10 TeV; otherwise, 
the BBN is threatened by the gravitino decay into the electronic and hadronic showers. 
Although we do not pursue this possibility, it is interesting to work in this directional 


3 Gravitino dark matter in 5D SUGRA 

3.1 4D effective Lagrangian and matter contents 

The soft SUSY-breaking terms are sensitive to the ultraviolet completion of the SM. As a 
concrete model, we consider the 5D SUGRA on S 1 / Z 2 and the flat 5D background metric, 

ds 2 = ri llv dx^dx v — dy 2 , 

where with /x = 0,1,2,3 and y denote the 4D spacetime and fifth coordinates, respectively. 
rjnv =diag(l, —1, —1, —1) and the fundamental region of the orbifold is chosen as 0 < y < L in 
which y = 0,L correspond to the fixed points. The S 1 /Z 2 orbifold restricts all fields f(x,y) to 
two classes of them such as Y 2 -even and -odd fields, satisfying the following Z 2 transformations 
f{x, —y) = f(x,y ) and f(x, —y) = —f(x,y), respectively. Only Z 2 -even fields have zero modes 
which can appear in the low-energy effective theory. 

First of all, we list the relevant matter contents of 5D SUGRA. From the structure of the 
orbifold, 5D SUSY is broken into the 4D J\f — l SUSY. Correspondingly, 5D vector multiplets 
V 1 and hypermultiplets are decomposed into 4D vector multiplets V 1 and three types of 
chiral multiplets Z 1 , <f> Q and that is, V 1 = {V 1 , S 7 } with / = 1,2, ,n v and H a = 
{<h a , } with a — 1, 2, • • • , Tin + nc where nc is the number of compensator hypermultiplets 

and in this paper, it is chosen as nc = 1, for simplicity. In addition to the usual Z 2 -even vector 
multiplets V 1 involving the vector multiplets in the standard model, we consider U(l)i> Z 2 -odd 
vector multiplets V 1 with I' — 1, 2, ■ • • , riy. The zero modes of Z 2 -even chiral multiplets E J are 
called as the moduli chiral multiplets T 1 ' whose linear combinatiorQ plays a role of the inflaton 
held as pointed out in Ref. |18j . In what follows, we define the zero mode of chiral multiplets < h Q , 
as Q a which involve the quark chiral multiplets {Qi,Ui,T>i), lepton chiral multiplets (£j,£j, Ni) 
with the number of generations i = 1,2,3, Higgs chiral multiplets {'H u ,'Hd), SUSY-breaking 
chiral multiplet A" and the stabilizer multiplets Hy. These multiplets have representations of 
the standard model gauge groups and extra U(l)p gauge groups whose gauge fields A : M , A : M 
in vector multiplets V 1 and V 1 ', respectively. It is then assigned 1/(1)// charges to these 
hypermultiplets H a . Here it is assumed that the visible sector consists of the MSSM plus 
right-handed (s)neutrinos and the same number of stabilizer hypermultiplets as that of moduli 
multiplets in order to generate the moduli and inflaton potential as can be shown later. 

3 We will discuss it in the separate work. 

4 In the case riy = 1, the radion multiplet corresponds to the single modulus T 1 =1 . 
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Next, we show the effective action obtained from the 5D conformal supergravity action for 
vector and hypermultiplets which is an off-shell description of 5D SUGRA |23, 24]. The struc¬ 
ture of the Kahler potential in 5D SUGRA on S 1 jZ<i can be characterized by the cubic poly¬ 
nomial of vector multiplets, the so-called norm function, Af(M) = YITjM j M k 
with real coefficients C^j^k for I,J,K = 1,2, ••• ,ny, and the U{l)r charges of hypermulti¬ 
plets. After the off-shell dimensional reduction discussed in Refs. [25, 26] [27] based on the 4D 
Af = 1 superspace [28] [29] JH the 4D effective Lagrangian is given by 


U e ff — — ~ 


d 2 0 fa(X, T)tr(W a W a ) + h.c. 


+ / d 4 0 |0| 2 O eff (|g| 2 ,ReT) 


d 2 0 0 3 W(Q,T) + h.c. 


( 4 ) 


where 0 is the compensator multiplet, W“ is the field strength supermultiplet for a massless 
4D vector multiplets V a with a = f/(l)y, SU(2)l, SU(3)c originating from the 5D 4 2 -even 
multiplets V a , Q a are the 4D chiral multiplets, X is the 4D chiral multiplet which induces the 
SUSY breaking, and T 1 are the moduli chiral multiplets. 

Then, the gauge kinetic functions f a (X,T ) in Eq. (J4|) are supposed as 


fly 

f.(X,T)=e x X + J2g, T 1 ', (5) 

I '=1 

where and f° x are real constants determined by the real coefficients Cpj^k hi the norm 
function and the gauge kinetic functions at the orbifold fixed point y = 0, respectively. Since 
the gauge kinetic functions at the orbifold fixed points depend on the dynamics of the SUSY- 
breaking sector, we comment on the reason why we take the above ansatz later. 

On other hand, the effective Kahler potential in Eq. (J4]) is given by 


n eS (\Q\ 2 ,ReT) = Af 1/3 (ReT) 


-S + 2J2 Y(c a ■ r)\Q a \ 2 + ^ o^(ReT)|g a | 2 |g^| 2 + o(|g| 6 ) 

a a,(3 

( 6 ) 


without the Kahler potential at the orbifold fixed points y — 0, L, where Af(ReT) is the norm 
function, Y(z ) = (1 — e~ 2Kez )/2Rez stands for the kinetic terms of Q a which have appeared 
after solving their equation of motion in the fifth direction and c : a denote the U(l)j/ charges of 
Q a . The four-point couplings 0^ are defined as 


q( 4 ) • Va- 1 • cp){Y((c a + c p ) • T) - Y(c a • T)Y(c p • T)} Y((c a + cy) • T) 

a 'P ~ (c Q -ReT)( C/ 3 -ReT) 3 

V j(X) = 


( 7 ) 


5 The more general 5D action, including i? 2 -odd fields, is discussed in Refs. (30] ]3l] 
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where V 1 j(X) is the operator to project the moduli multiplets out the radion multiplet. The 
notable feature there is that the flavor structure of matter fields is characterized by the U{l)p 
charges of them in the Kahler potential ([6]). By contrast, the superpotential can be allowed 
only at the orbifold fixed points where the SUSY is reduced to the 4D J\T = 1. Therefore, we 
consider the superpotential including the Yukawa couplings and //-term in the MSSM, moduli 
potential at y — 0, and the moduli potential at y = L, respectively. The explicit form of the 
superpotential in Eq. (J4]) is shown later. 


3.2 Gravitino dark matter in 5D SUGRA 


In this section, we show the realization of mass hierarchies between the gravitino and other 
sparticles in the framework of 5D SUGRA. As shown in Eq. (JBj) in Sec. 13.11 the bulk Kahler 
potential is rewritten as 

Abuik = -lnAA(ReT) + ^ Z Qa (ReT)|Q a | 2 + Z v (ReT)|Y| 2 + 0(|Q| 4 ), (8) 


where the Kahler metric K xx for the SUSY-breaking held X is given by 


K X x = Z x (ReT) 


_ g—2cx-ReT 

cx ■ ReT 




cx -ReT ’ 

^^|exp(2|cx-ReT|) 


c x ■ ReT > 0, 
cx ■ ReT < 0, 


(9) 


where K xx depends on the U(l)p charges of the held X for the Z 2 -odd vector multiplets V 1 ' 
and the VEVs of the moduli T 1 , except for the case of the vanishing U(l)i> charges. For the 
mild large volume of the hfth dimension, L ~ A^ 1//2 ((ReT)) 1 and positive 17(1)// charges, 

the VEV of the Kahler metric is smaller than 0(1), that is, (K xx ) -C 1, which is important to 
obtain so that the light gravitino can be lower than the other sparticles. 

The soft SUSY-breaking masses for the scalar components of Q a are given by the four-point 
couplings in Eq. (J7J). For typical U(l)p charges of Q a to realize the realistic Yukawa 

couplings, the soft SUSY-breaking masses are larger than the gravitino mass as shown later. 
Furthermore, the gauge kinetic functions in Eq. (J5J) lead to the following gaugino masses at the 
compactihcation scale by employing the formula (l3j) , 


M a 


F x 
9 a 



I '=1 



( 10 ) 


When the compactihcation scale is close to the GUT scale, we obtain the gaugino masses at 
the EW scale after solving the one-loop RG equations from the GUT scale to the EW scale, 


Mi (Mew) — 0.4M!(Mgut) ; M 2 ( M ew) — 0.8M 2 (Mgut) ; ^7 3 (Mew) — 2.9M 3 (Mgut)- (H) 


Then, the gravitino LSP occurs if these gaugino masses at the EW scale are larger than the 
gravitino, as pointed out in Ref. [21 j . In the case of 5D SUGRA, such situations can be realized 
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by properly choosing the parameters and at the same time, the Higgsino mass should be 

larger than the gravitino mass. Thus, one can consider the gravitino dark matter in the gravity 
mediated SUSY-breaking scenario without changing the VEVs of the F-terms as discussed in 
Sec. EJ In order to estimate thermal and nonthermal abundances of gravitino via the moduli 
and/or inflaton decay, we focus on the specific model which realizes the successful moduli 
inflation as well as the moduli stabilization [18] in the next section 13.31 

The mild large volume also reduces the contribution from the Kahler potential at the orbifold 
fixed points y = 0,L to be small compared with the bulk Kahler potential ()8|) . Since their 
boundary terms are described by 

^boundary = A^ 1 / 3 (a^(|X| 2 ) + (e~ C ^ T+ ^\X\ 2 ) + ■■■), (12) 

the overall factor Af~ 1//3 suppress these contributions. The one-loop corrections to the moduli 
Kahler potential [32] are also suppressed by the mild large volume of the fifth dimension. 

By contrast, in the case of negative U( 1)// charges, the VEV of the Kahler metric is bigger 
than 0(1), that is, (K X x) ^ 1- From the mass formula of the gravitino and sparticles given 
by Eqs. ([2]) and (Ffl) , one can expect that the sparticles are lighter than the gravitino without 
changing the F-terrn of the SUSY-breaking field. Thus, it is possible to solve the gravitino and 
fine-tuning problems at the same time. 

3.3 Moduli stabilization 

Following the discussion about the small-field inflation in Ref. [T8] , we choose the norm function 
as 

AZ’(ReT) = (ReT 1 )(ReT 2 )(ReT 3 ), (13) 

which leads to the diagonal moduli Kahler metric. Because it seems to be difficult to obtain 
the realistic masses and mixings of quarks and leptons in the case of two moduli as shown in 
Sec. 14.11 we restrict ourselves to the case of three moduli T 1 =1,2,3 in what follows. In order to 
generate the moduli potential, we introduce the same number of stabilizer chiral multiplets Hi 
as that of moduli chiral multiplets as stated in Sec. 13.11 The effective Kahler potential, except 
for the SUSY-breaking field A" and other matters in the MSSM are 

3 

K = —lnJ\f(ReT) + ^(ReT J ' =i )|Hi| 2 , ( 14 ) 

i— 1 

where it is then assumed that the stabilizer fields Hi have only the U{l)r=i charge with i = 
1, 2, 3, for simplicity. In addition to it, the relevant superpotential for the moduli inflation and 
stabilization is 

iu,od = e YU” - E 1 

i= 1 i=l 

= E (A - as) 

i= 1 
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where Jq\ j are constants at the orbifolcl fixed points y = 0, L and the exponential factor e c p T 

comes from the profile of the wavefunction of the stabilizer fields in the fifth direction, ^ = 
(o') 

e i 7 Here we assume that these tadpole terms are dominant in the superpotential and 

the other terms are negligible due to some symmetries or dynamics^ In the following, we omit 
the subscripts of the stabilizer fields at the fixed point y = 0, that is Hi = H 

In fact, from the 4D scalar potential (CQ) given by the Kahler and superpotential (1X41) . (fT51) . 
the expectation values of the moduli T 1 and the stabilizer fields Hi are found as [19] 

cf{T e ) = In 4 (H t ) = 0, (16) 

Jo 

which are determined by the stabilization conditions, (Dj/W) = ( DiW) = (W) = 0, at which 
the supersymmetric Minkowski minimum can be realized, (V) = 0. Their supersymmetric 
masses of moduli and stabilizer fields are estimated as 


m 


2 

Pi 


e {K \W ri ) 2 

(K Iip )(K a y 


(17) 


where ( Wn) = — c^./}^e _c n T and Wij = didjW. Now there are no mixing terms between the 
moduli and stabilizer fields in the mass matrices due to the diagonal Kahler metric of them. 
From the exponential behaviors of supersymmetric masses (ITT]) , the mass scales of moduli and 
stabilizer fields are controlled by the sizes of U(l)p=i charge and constants Jq L . 

So far, the SUSY is not broken in the superpotential (1T5]h For the SUSY-breaking sector, 
we consider the O’Raifeartaigh model [3X] which is simplified as the following Kahler and 
superpotential of the SUSY-breaking field X after integrating out the heavy modes, 


K = Z x (ReT\ReT 2 )\X\ 2 - ^-\X\\ W = w + uX , (18) 

A- 

where w, v are the real parameters and the SUSY-breaking field X has no U(l )3 charge, for 
simplicity. The Kahler potential receives the loop corrections from the heavy modes, whose 
mass scale is A [35] . 

In general, the true vacuum of the moduli and stabilizer fields are deviated from the su¬ 
persymmetric one due to the SUSY-breaking effects and then the moduli and stabilizer fields 
obtain their F-terms at the true vacuum. Since their F-terms would change the cosmological 
history of the universe through the moduli decay into the gravitinos, it is important to evaluate 
their F-terms at the true vacuum. For that reason, we adopt the perturbation method, known 
as the reference point method [36] to search for the true vacua of all the fields. 

First, as the reference points for the moduli and stabilizer fields, we take them as given in 
Eq. (1T6|) satisfying as 


D H W |ref = W Hi + K H W 


0, F> r /'IU| re f = K tI 'W , 


(19) 


6 A similar moduli stabilization was proposed in Ref. j33i in the case of nc = 2. 




and for the SUSY-breaking field A", its reference point is taken as that satisfying the following 
stabilization condition: 


e~ K V x \ rei = d x (J2 KTl ' fI '\D t i'W\ 2 + K XX \D X W\ 2 - 3\W\ 2 ) 
v 

~ 3 W X W + d x {K xx )\W x \ 2 + K xx W x K xx W - 3 W X W 
4|hUx| 2 - 

~ ———A + WvW = 0 
“A 2 (Z *) 2 ’ 

in the limit w 1, where V x = d x V. Thus, we obtain 

A 2 (Z.v ) 2 ( W \ A 2 (Zx) 2 w 

X ^ = —— (wyj"— — 


( 20 ) 


( 21 ) 


Next, we expand these fields as 0 —* 0| re f + 50, 0 = T v , Hi , X with V, i — 1,2,3 and evaluate 
their perturbations from the reference points given by Eqs. (TT6l) and (T2TT) under the following 
conditions: 


V = U| ref + V] | ref 50^ + VjUM 1 + U 7J | ref 50 7 50 J + Vj j\ re { 6(f) 1 5(j) J + U 7j -| ref 50'50 J + O(50 3 ), 


v^M 1 + VrUsp 


> 


V/j| ref (50 / 50 J + U J j| ref 50 / 50 J + V- I j\ rei 5(j) I S(t) J 


( 22 ) 


where Vj = djV and Vu = didjV are the derivatives for the relevant fields 0, and then 0| re f + 50 
are considered as the vacua of relevant fields. Note that these perturbations are valid even if 
the SUSY-breaking scale is smaller than the scale of supersymmetric masses given by Eq. mm- 
As a result, the deviations of the fields from the reference points (12TT) are 


A fj ~_ 

1 ~ 2ReT I 'W T p H . ’ 



5X ~ 


4 w 2 


5wW x , 


(23) 


and the E-terms and squared masses of moduli, stabilizer, and SUSY-breaking fields are roughly 
estimated as 


Vk, 


j'l'j'l' 


f t ' ~o 


w 


m, 


1 j 1 !' 

e K W 2 Tl , Hi 


— m 2 H . — K K 




(!' = i), rr>x ^ 


w 


m 


Hi 


, \/~KxxF X — 




e K Aw 2 


Kxx^zy 


(24) 


at the vacuum, 0 = 0| re f + <50. The mass squares of real and imaginary parts of moduli, 
stabilizer, and SUSY-breaking fields are the same as each other and here and in what follows, 
they are denoted as , m 2 Hi and m 2 x , respectively. The details of these derivations are 
summarized in Appendix. [A] The mass differences between m T i> =i and m h, are the order of the 
gravitino mass. It is remarkable that the fields, except for the SUSY-breaking held X, have 
almost vanishing E-terms due to their large supersymmetric masses. 
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3.4 Moduli inflation 

In this section, we briefly review the inflation mechanism in which the inflaton is identified 
as one of the real parts of the moduli. Although in Ref. [18], both the small- and large-field 
inflation are discussed in the light of recent Planck results, in this paper, we restrict ourselves 
to the small-field inflation, for simplicity^ The inflaton potential is generated by the Kahler 
and superpotential of the pair (T 3 ,ii 3 ) in Eqs. (IT-fll and (fT5]h 

K = -hiA/'+Ztf 3 (ReT 3 )|H 3 | 2 , 

bhinf = (4 3) - 4 3) e- c 3 3)t3 ) H :il (25) 

where 

1 _ —2cg 3) ReT 3 

z - (ReT3) = -^5^-- (26) 

and the effective scalar potential is obtained from Eq. (P) with the above Kahler and superpo¬ 
tential (1251) . 


R nf = e K K H ^\W H3 \ 2 ~ 


I J, 


(3) 


j(V4 3 >T3|2 


(ReT 1 )(ReT 2 )(l-e-4 3) ^) 


T 3 > 


(27) 


where ReT 3 is identified as the inflaton. Here, it is supposed that the other moduli T 1 ', 
stabilizer fields Hi with T, i = 1,2 are heavier than the pair (T 3 , H s ) and fixed at their minima. 
This is because the minima of them are fixed by their own superpotential in Eq. (1T51) . and they 
can be decoupled from the pair (T 3 , H 3 ) by choosing the parameters in the superpotential (1T51) , 


j W — j 

Jr\ — d ( 


( 2 ) 




(28) 


and the nonvanishing (1) 1 , 2,3 charges of i/ 1 , 2 , 



1 

50’ 


(29) 


whereas the constants are chosen to be small compared with as shown later. Fur¬ 

thermore, in the following analysis, we omit the fluctuation of H 3 and A", because their minima 
are fixed around the origin by the Hubble-induced mass during the inflation. Irn T 3 is also fixed 
at the origin during and after the inflation. They can be checked that the fluctuations of these 
fields are negligible to the inflaton dynamics as explicitly shown in Appendix. iBl 

When Re T 3 is identified as the inflaton, the effective scalar potential (f27|) is similar to the 
one in the Starobinski model [37] and is drawn in Fig. |T] with the parameters given by Eqs. (1251) . 
(129|) and (j36|) . From Fig.P the inflaton, ReT 3 can roll its potential slowly down to its minimum 

7 The extension to the large-field inflation is straightforward. 
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from the large value of ReT 3 . In order to evaluate the cosmological observables for the cosmic 
microwave background (CMB) observed by Planck, we define the slow-roll parameters for the 
inflaton, a = ReT 3 , 


e = 


Mp\ 

2 v;i, 


V = 


V <J V (7 Idnf 

VW 


(30) 


where V CT is the Kahler covariant derivative for the held a. With these slow-roll approximations, 
the power spectrum of the scalar curvature perturbation, its spectral index, and tensor-to-scalar 
ratio can be expressed as 


m) 


1 hjnf 
247r 2 e ’ 
dlnp^k) 

dink 


~ 1 — 6e + 2 V , 


r = 16e. 


(31) 


The recent data reported by the Planck Collaboration shows the almost scale invariant spectrum 
and the upper limit of r [22j , 

Pt(k) ~ 2.20 ±0.10 x 10" 9 , n s = 0.9655 ± 0.0062, r < 0.11, (32) 


at the scale k * 


0.05 [Mpc 1 ]. The inflaton dynamics is obeyed by its equation of motion, 


// 

a 



3 a + 6a 2 4r ) + ^ 


.V 

V 


a 


(33) 


where ' denotes the d/dN by employing the number of e-foldings N rather than time; 

. . M d dN d d , d . 

“(«) = 17 = -J7-T ,7 = H-7Tr. ( 34 ) 


where a(t) is the scale factor of 4D spacetime. The metric g aa is connected to the Kahler metric 
K t 3 f3 such that | g aa dodo = K T 3f3dT 3 dT 3 and Y a a(J = —1/a is the Christoffel symbol. As a 
result, it is found that the power spectrum of scalar curvature perturbation, its spectral index 
and tensor-to-scalar ratio are consistent with the current cosmological data, 


P ( (k) ~ 2.2 x 1(T 9 , n s ~ 0.96, r ~ 10" 5 , 


(35) 


with the enough e-foldings N ~ 58 and then the parameters are chosen as, 

t( 3 ) _ 1 t(3) _ 3 (3) _ J_ 

0 4000’ L 4000’ 3 10' 


(36) 


The running of the scalar spectral index is negligible, relative to the current observational 
sensitivity. 
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Unf 



Figure 1: Inflaton potential Unf on the Im T 1 = 0 hypersurface. 

When the numerical values of the parameters are chosen as those in Eqs. (128|) . (1290 and 
(T36|) . the moduli VEVs become 

(Re(T 1 ), Re(T 2 ), Re(T 3 )) ~ (110, 110, 11), (37) 

in the unit Mp\ = 1. According to these moduli VEVs, the typical Kaluza-Klein mass scale is 
found as 

Mc = i~ - 2X X 1Ql6GeV ’ ( 38 ) 

which is close to the GUT scale due to the mild large volume of the fifth dimension, ( JV 1//2 ) ~ 
364. The mass of moduli T 1 and stabilizer fields are also given by 

m T i ~ m T 2 ~ m Hl — m H2 ~ 4.8 x 10 15 GeV, m T 3 ~ m h 3 ~4x 10 12 GeV, (39) 

and their F-terms become 

F rl ~ F t2 ~ F Hl ~F ft ~lx 10” 42 , F t3 ~ F Hs ~ 1.6 x 10 -36 , (40) 

in the unit Mpi = 1. 

So far, we have specified the parameters relevant for the moduli and stabilizer fields. The 
parameters in the Kahler and superpotential (fT8]( for the SUSY-breaking sector are considered 
as 

z/ ~ -1.567 • • • x 10“ 14 , w = —6 x 10' 14 , A = 10" 4 , = ^, cj } = c^ = 0, (41) 

where v is proper chosen as realizing the Minkowski minimum. Then, the mass of the gravitino 
and the mass and F-term of X are obtained as 

F x 

m 3 /2 — 395 GeV, nix — 6 x 10 8 GeV, —— ~ 4541 GeV, (42) 

Mpi 

which implies that the gravitino mass is suppressed by the Kahler metric of the SUSY-breaking 
field, K X x — 1 /(cx ■ ReF) ~ 0.023 as discussed in Sec. 13.21 and the concrete sparticle spectra 
are shown in Sec. 14.31 
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3.5 Moduli-induced gravitino problem and reheating temperature 

As mentioned in Sec. 13.31 the moduli and stabilizer fields are so heavy that they decay into the 
particles in the MSSM before the BBN. However, even if they are much heavier than 0(100 
TeV), it has to be taken into account of the cosmological problem, e.g., the moduli-induced 
gravitino problem [38], 139]. 

The moduli decay width into the gravitino pair can be evaluated by the couplings between 
moduli and gravitinos in the unitary gauge, 

£ 3/2 = - e K / 2 W*^a^ v - e^W^a^, (43) 


where p p denotes the gravitino in two-component formalism and the relevant covariant deriva¬ 
tives of the gravitino are D p p a = d p p a + - l (Kjd p (p J — Kjd p (p J ). After carrying out a field- 
dependent chiral transformation, 





(44) 


the Lagrangian (1451) is simplified as 


^fiispcr 

C 3/2 = ^ vp<J ^pP v d p p a + — —(Gjd p (j) J - - e G/2 {p p o yiv p v + $^$„), ( 45 ) 

where G = K + ln|iy| 2 and Gj = djG. When we expand the moduli T 1 ' around the vacuum 
given by employing the reference point method (fT6l) and (1251) . the Lagrangian (1431) reduces to 

e /4i/ptr _ _ 

£ 3/2 =-^yrr.Vyl', + —({G TJI )d f ST JI - (G T j,)d p ST 

- - 1( e G ' 2 )((G TJ ,)ST J ' + {G-j.)ST r )$ |1 |r,7‘1®'., (46) 

in the four-component formalism of the gravitino X I J p . As shown in Eq. (1391) . the moduli and 
stabilizer fields, except for the pair (T 3 ,H 3 ), are decoupled from the inflaton dynamics due to 
their heavy masses. Therefore, their decays can be neglected and do not give the sizable effects 
in the thermal history of the universe. In this respect, we focus on the decay processes of T 3 , 
H 3 , and SUSY-breaking field A". 


3.5.1 The inflaton decay 


First, we concentrate on the inflaton decay into gravitino pair. Since the gravitino wavefunction 
is described in terms of helicity ±1/2 components of the gravitino at a high-energy limit by the 
equivalence theorem, the inflaton decay width into the gravitino pair is estimated as 


T(a 3 —>- Ts/od^) — 


1 

2887i(K T 3f3) 


D t3 W 

w 


~ 1.6 x 1(T 18 GeV, 


m b T 3 1 rn T 3rnlj 2 

— 288ti (K T 3f3N) 1 

(47) 
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by employing the E-term of modulus (1241) and numerical values of the mass, E-term of inflaton, 
and gravitino mass given by Eqs. (T39|) and (140|) . Here, the reduced Planck mass has been explic¬ 
itly written. When the inflaton has the sizable E-term at the vacuum, the enhancement factor 
m^/ 2 2 , as the longitudinal mode of the gravitino, induces the significant amount of gravitinos 
which would threaten to destroy the success of BBN. However, in our moduli inflation, this 
direct decay is so suppressed due to the almost vanishing E-ternr of the inflaton. Therefore, 
the dominant decay process of inflaton comes from the interactions with the gauge bosons, 


^Tgg ~ 


4(<?c 


_ T^a r^afiu 

\2 r 


J' 


-fj'STk 


T^a rpafiu 


-Q,ST}' F‘F‘ 


(48) 


where = ReT J/ , Tf = ImT^. Now, the gauge kinetic functions f a (X,T ) are considered 
as in Eq. (J5]). In general, ^ a x X in the gauge kinetic function ([5]) could appear, because the 
E-symmetry is explicitly broken by the constant superpotential w in Eq. (fT 8 |) . 

The inflaton decay width into the gauge bosons are 


3 

£ 

a= 1 


IV 


->■ g (a) g {a) ] 


n g / a y 

1287T \ y/2K T 3f3 I 


m 


T 3 


Ml\ 


~ 3.95 GeV, 


(49) 


with the numerical values of mass and VEV of modulus fj37jl . (139|) . where Nq are the number of 
the gauge bosons for the gauge groups U(l)y, SU( 2 ) L , SU (3)c and the nonvanishing coefficients 
in the gauge kinetic function are chosen as £3 = £3 = = 0.22 to realize the gauge coupling 

unification at the GUT scale Mqut — 2 x 10 16 GeV. Although there are the other decay 
processes via the inflaton decay into the gauginos A a given by the interactions, 


Ctxx = —^Re/ a (A°cffiZyv + (H.c.)) + ^Im/ a D # 1 (AW“) 

a a 

+ £ + (H.c.)Y (50) 

where D M is the covariant derivative for gaugino, such decay channels are suppressed by the 
small masses of gauginos and almost vanishing E-term of inflaton such as 

rV -> A“A«) ~ y ~ 1,5 x 10" 21 GeV. (51) 

v ^ 16 tt Ml, 

a= 1 

with m \3 ~ 1.5 TeV and the derivative of E-term for the inflaton, 

(l^) = {wX' 2 + K™D T ,WD^W)) ~ O (J) . (52) 

The decays from the inflaton into sfermions are also suppressed because of the factor, 
^sfermion/wp, if the masses of sfermions are of 0(1 TeV). Other decays from the inflaton into 
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the fermion pairs and quark-quark-gluon are negligible due to their small masses and phase 
factors, respectively, as pointed out in Ref. [39] . The /i-term does not give the sizable effects 
for the inflaton decay process, because we consider the tiny /i-term (~ 500 GeV) in the light 
of naturalness as shown in Sec. 14.11 Finally, we comment on a single gravitino production via 
the inflaton decay into the modulino and gravitino. Since the mixing terms between T 3 and 
H 3 in the mass squared matrices are controlled by the SUSY-breaking scale, i.e., the gravitino 
mass, the mass difference between the inflaton and modulino as its superpartner is of the order 
of the gravitino mass. Therefore, the inflaton decay width into the modulino a 3 and gravitino 
is suppressed by the phase factor m 3 /2/m T 3 , 

r(ff 3 _> rt 3/2) * Y (^) 2 »s /2 - 7.2 x 10-GeV, (53) 

with 777 , 3/2 = 395 GeV, m T 3 = 4 x 10 12 GeV given by Eqs. (l39]i and fl42]h The inflaton decay 
into the SUSY-breaking field X is also suppressed, because there is no tree-level interaction 
between X and T 3 due to the vanishing { 7 ( 1)3 charge of X. As a result, the branching ratios 
of the moduli decaying into the gravitino(s) are summarized as 

3 

T^i = T(cx 3 all) ~ J^T(a 3 ->• g^g {a) ) ~ 3.95 GeV, 

a —1 

Br( ff 3 -5 * 3/2 * 3/2 ) c F(g3 c 1.4 x 10-“, 

all 

Br(cr 3 ->■ d 3 T 3/2 ) ~ ^ 3/2 ^ - 1.8 x 10 -22 , (54) 

all 

and then the reheating temperature is roughly estimated by equaling the expansion rate of the 
Universe to the total decay width of inflaton, 

/ 2 \-V 4 

F& = H r &T r = VraJiMpi ~ 1.38 x 10 9 GeV, (55) 

where Hr = H(Tr ) and g*(Tn) = 915/4 is the effective degrees of freedom of the radiation in 
the MSSM at the reheating. The gravitino yield I 3/2 via the inflaton decay is suppressed due 
to the tiny branching ratio of the inflaton decay into the gravitino(s), 

Y 3/2 = ^ ~ Br(cr 3 -> cr 3 'F 3/2 )-p—~ 3.8 x 10" 24 , (56) 

s Ahits 

with m 3 /2 = 395 GeV, Tr = 1.38 x 10 9 GeV s = 4p/3T, where n 3 / 2 , s, and p are the number 
density of the gravitino, entropy, and energy density of the Universe, respectively. Now it is 
supposed that the coherent oscillation of the inflaton field dominates the energy density of the 
Universe after the inflation and there is no entropy production after the inflation as shown 
later. 
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It is remarkable that the supersymmetric moduli stabilization is important to suppress the 
direct decays from the inflaton into the gravitino(s) which give the solution to the cosmological 
moduli problem, especially the moduli-induced gravitino problem. The other gravitino produc¬ 
tion from the stabilizer fields, the SUSY-breaking held, and the thermal bath can be estimated 
in the next section. 


3.5.2 The decay of stabilizer and SUSY-breaking fields 

The stabilizer field H 3 is stabilized at the origin during the inflation and after that, Re H 3 
oscillates around its vacuum (123|) deviated from the supersymmetric one (1161b On the other 
hand, Im H 3 and Im X evolve to the origin during inflation and do not oscillate after the inflation 
as shown in Appendices. [A]and|B] Similarly, the SUSY-breaking field Re A" oscillates around 
its vacuum after the inflation. From the analyses in Appendices. |A] and [Bl the amplitudes of 
both fields are found as 


A h 3 - 


m 3/2 

rriH 3 ’ 



(57) 


with h 3 = Re H 3 and x = Re A". By comparing their masses given in Eq. (l39]i with the reheating 
temperature (l55lh the coherent oscillations of both fields h 3 and x start before the reheating 
process. When H 3 does not couple to the fields in the MSSM, the dominant decay process is 


rjs = r (h 3 t 3 /2^3/ 2 ) ^ 


288t tK, 


1 h 3 h 3 
-0.02 GeV, 


Dh 3 W 

W 


m 


h 3 


m 


h 3 


ml /2 Ml L 2887 t(K H3S3 M) M& 


(58) 


which implies the decay time of h 3 is smaller than the time of the coherent oscillation of h 3 
and reheating, that is, H^ c > H(T R ) > H^ c , with H^ c - m h3 and H^ c - T(h 3 -> ^ 3/2 ^z/ 2 )- 
Here and in what follows, H R , H® sc , and Hf ec refer to the Hubble parameters at the time of 
reheating, beginning of oscillation of relevant fields <h, and decay of <h. The scale factors of 4D 
spacetime a R , a® sc , and af ec are also defined in the same way as the Hubble parameters, H R , 
H® sc , and Hf ec . The energy density of coherent oscillation h 3 is 

3 , (59) 

" ^ \ ^OSC / 


where a^ c stands for the scale factor at the time when h 3 begins to oscillate and ph 3 is converted 
into the gravitino yield hereafter, 



2 p h3 _ 1 m l/2 T R 
•m Hi s 4 rriff 3 


- 8.2 x 10~ 25 , 


(60) 


with m 3 /2 = 395 GeV, T R = 1.38 x 10 9 GeV, and m H3 = 4 x 10 12 GeV. Here we employed 
that the entropy production from h 3 can be neglected. In our model, the following inequality 
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is satisfied due to the tiny mass of the gravitino and then h 3 does not dominate the Universe 
and release the significant entropy, 


1 > 


Ph 3 

Pr 


T=T f 13 

dec 


Ph 3 


( M 

P 

end 

luJ 


m 3/2^Pl / Tr \ 

2 Vjnf \ T h 3 J 

\ r dec / 


(61) 


where ph 3 , Pr are the energy densities of h 3 and radiation, respectively. p| en d = Unf — 0(10 13 ) 
denotes the energy density at the end of inflation and T^ c is the decay temperature of /i 3 , 


_ 

1 dec 



-1/4 

yrJgMp, ~ 8.6 x 10 7 GeV. 


(62) 


Furthermore, the SUSY-breaking field also produces the gravitinos through the following 
dominant decay channel: 


T(x ->■ ^ 372 ^ 3 / 2 ) - 


1 

/ D X W \ 

2 m x 1 

V 

2 m x 

288tt(K xx ) 

\ W / 

2887 t(K xx ) 

w 



(63) 


With the parameters (Hill , the VEVs of moduli (j37jl . and mass of 1113/2 and X (l42lh the total 
decay width of A" then becomes 


rj, = r(i -> ® 3 / 2 ® 3 / 2 ) - 3.7 x 10 s GeV. 


(64) 


Therefore, the decay time of x is smaller than that of reheating, that is, H* sc > «S* » H(T r ), 
with H£ sc ~ rrix and H^ ec ~ T(x —> ^ 3 / 2 v I , 3 / 2 ) - The energy density of the coherent oscillation 
x is converted into that of the gravitino as 


Px 





1 

2 




(65) 


at the time of reheating, where the gravitino is relativistic at the time of production. By 
employing the scale factors, 


Rr 


'Ven 


TUrS 


the gravitino yield is 


vx _ 
X 3/2 — 


-2/3 


y/6m x 

771r3 


- 2/3 


“dec 
’ rr’3 

a* 


6 ( r all" i A ') 5 


-1/3 


m 


( 66 ) 


2 p x ^ 3 T r / m 3 / 2 
m x s 2 nix V m x 


16/3 


all 


;\ 2/3 

!■ J ~ 2 x 10 


3/2 


-32 


(67) 


with the numerical values given by Eqs. (J42J), (J54]), ((55]), and ((64]). It is found that the gravitino 
production via x decay is suppressed by the tiny mass of the gravitino and it is not the dominant 
source for the relic abundance of the gravitino. The entropy production from x can be also 
neglected in the same way as that of h 3 . As pointed out in Ref. [3D|, under 7713/2 m x -C 
m T 3 < A, the gravitino production is signihcantly relaxed and this condition is satisfied in our 
model. Note that when A is smaller than the inflaton mass, we have to take account of the 
inflaton decay into the fields in the hidden sector. 
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4 Gravitino dark matter and the Higgs boson mass 


4.1 Yukawa couplings and naturalness 

Before estimating the relic abundance of the gravitino, we specify the Yukawa couplings and 
/i-terrn in the superpotential which can be only introduced at the orbifold fixed points y = 0, L, 
where the SUSY is reduced to 4D AT = 1. As stated in Sec. 13.11 we consider the Yukawa 
interactions in the MSSM at the orbifold fixed point y = 0, 


UAukawa — KjQiHuUj + A fjQi'Hd'Dj + A i j £ i 'H d £j + X'^C.H.j.Xj. (68) 

where A7 d,e,n are the holomorphic Yukawa coupling constants and are supposed to be of 0( 1). 
After the canonical normalization of fields in the MSSM, the physical Yukawa couplings are 
expressed as 


Vij = 


\U 

\j 


VWeJnM)’ 


4 


\d 

*ij 


^/(Y Q Y Hd Y Vj y 


Va = 




vlW 


Vij = 


\n 

\j 


(69) 


where 

Y a = 2Y 1 / 3 (ReT) |V(c a • T) + 0^(ReT)|Y| 2 + 0(|Y| 4 )} (70) 

The function Y(z) is always positive, and approximated as 


Y(z) 


_ g—2R ez 

2R ez 




2Rez ’ 

«exp(2|Re^|). 


Rez > 0 
Rez < 0. 


(71) 


In the 5D viewpoint, the wavefunctions of fields are localized toward y — 0 (y — L) in the 
case that c a ■ (ReT) is positive (negative). As can be seen in Eq. (ITT]) . y^- d,e,n are of 0(1) or 
exponentially small when all the relevant fields are localized toward y = 0 or y = L, respectively. 

Therefore, we expect that the mass hierarchies of elementary particles and the extreme 
smallness of the neutrino masses can be realized even in the case of Dirac neutrinos. In fact, 
when we choose the (7(1)// charges and 0(1) values of the holomorphic Yukawa couplings 
yu,d,e,n - n ^pjgg Q] arK ] [2j the observed masses and mixing angles of quarks and leptons at 
the electroweak scale can be realized. Here, we employ the full one-loop RG equations of the 
MSSM from the GUT to the EW scale. It is remarkable that the flavor structure of soft SUSY- 
breaking terms is determined by the (7(1)/' charge assignment as can be seen in the Kahler 
potential dSJ). In fact, the soft SUSY-breaking terms at the GUT scale are determined by the 
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following formula (4T, 02] : 


M a = (F'd ,In (Re/“)>, 
m Q. = -{F'F J drdj\nY Ql ,), 

A% = {F‘d,ln(Y H „Y ei Y Ul )), 

a % = (F'S/in (y^y-a,^)), 

At j =(F I d I \n(Y Hd Y Cl Y £l )), 

A"j = (F’dM (Y n „Y Ct Y Ml )), (72) 

where indices I and J run over all the chiral multiplets. Then, the [7(1)/' charge assignment 
in Table [Q and the F-terrn of the SUSY-breaking field X given by Eq. (1421) give rise to the soft 
scalar masses and gluino mass in Table [3] By contrast, the H-terms are almost vanishing due 
to the tiny F-terms of moduli. Here and hereafter, we parametrize the ratios of gaugino masses 
at the GUT scale as 

_ Mi(Mgut) _ M 2 (Mgut) 

M 3 ( M G ut ) ’ M 3 ( M G ut ) ’ 

where Mi(Mgut), M 2 (Mgut), and M 3 (Mgut) are the bino, wino, and gluino masses at the 
GUT scale, Mqut — 2 x 10 16 GeV. The ratios of gaugino masses are controlled by the param¬ 
eters £y in the gauge kinetic function (J5]) without spoiling the gauge coupling unification due 
to the tiny VEV of X. 


4f = (0.1, 0.1, 1.1) 

c q~ 2 — (—0-1, —0.1, 0.8) 

4= 3 = (0.1, 0.4, 1) 

4'r 1 = ( 0 . 1 , o.l, 1 . 6 ) 
c cf = (-0-1,-0.1, 0) 

4'r 3 = (0.1, 0.5, 0) 

1 = 0 
riu 

Cy =2 = 0.1 
rtu 

c^= 3 = -0.9 

riu 

eft 1 = (0.1, 0.1, 0.6) 

c u~ 2 ~ (—0-1, —0.1, 0.3) 
4= 3 = (-0.2, 0.2, 1) 

4r l = (0.!, 0.2, 0.2) 

4r 2 = (-o.i,-o.i,o) 

4= 3 = (-0.2, 0, - 0.5) 

enf = 0 

4r 2 = 0 

4/f = -o.i 

c T '= x = (0.1, 0.1, 0.2) 
4= 2 = (-0.1,-0.1,0) 
4= 3 = (0.3, 0.2, - 0.5) 

= (0.1, o.l, o.l) 
c 7 y= 2 = (-0.3, -0.3, -0.3) 

4= 3 = (-0.7, -0.7, -0.7) 



Table 1: U{l)r flavor charges of the quarks, leptons, and Higgs. 


On the other hand, the /r-terrn can be generated by the following superpotential: 

3 

Wfi -term — E KiHiHuHd, (74) 

i =1 

where K t are the 0(1) dimensionless couplings, Hi are the stabilizer fields with F-charge 2, 
whereas Higgs chiral superfields do not have the F-charge. These cubic interactions do not 
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w>l 

\K\ 

1 

( 0.32 0.35 0.95 ^ 
0.22 0.42 0.33 
y 0.51 0.48 1.5 ) 


1 

( 0.45 0.5 0.59 ^ 

0.28 0.24 0.38 
^ 1.03 1.02 0.81 ) 


|A7| 

IASI 

1 

( 0.28 0.22 0.52 ^ 
0.4 1.15 0.31 

^ 0.8 1.02 1.05 ) 


1 

( 0.77 0.85 0.69 \ 
0.25 0.98 0.58 
^ 0.34 0.26 1.03 ) 



Tabic 2: 


0{ 1) values of the holomorphic Yukawa couplings \^- d,e,n 


in the superpotential fj68jl . 


affect the moduli stabilization as well as the moduli inflation due to the almost vanishing 
VEVs of the Higgs fields. Thus, it is possible to consider the VEVs of the stabilizer fields Hi 
as the origin of the fx- term. After the canonical normalization of the relevant fields, the fx -term 
at the GUT scale becomes 

11 ~ t! (ubvTU 

Especially, in the case of k 2 = = 0, the scale of the /i-terrn is chosen as TeV scale, 



H ~ 3.8 x 10- 3 ^^M P1 ~ 0(500 GeV), 
m Hl 


(76) 


where K\ = 2/3, m 3 / 2 = 395 GeV, m Hl ~ 4.8 x 10 15 GeV, and (Hi) ~ rnz/ 2 /m Hl are given by 
Eq. (123|) and the factor 3.8 x 10 -3 comes from the mild large volume of the fifth dimension and 
normalization factors for H \, H u , and H ( i . The EW symmetry breaking requires the following 
relation between the mass of the Z-boson, mz and soft SUSY-breaking masses of the up-type 
Higgs m Hu : 


~Y - - |fl(M EW )| 2 + 0 . 


(77) 


in the limit of large value of tan (3, where /i(M EW ) and m Hu (MEw) are the /i-term and mn u 
at the EW scale, respectively. The VEVs of up- and down-type Higgs fields are denoted by 
v u = v sin/3 and Vd = v cos/3 with v = 174 GeV. Thus, the observed Z-boson mass indicates 
|/r(M EW )| ~ ~ otherwise /i and mn u have to be properly tuned to obtain the 

EW vacuum. We adopt the measure of the degree of tuning the /i-terrn at the GUT scale as 


A, 


1 d In m 2 z 

2 d In \fx\ ’ 


(78) 


and then 100 x % represents the degree of tuning to obtain the Z-boson mass mz = 91.2 
GeV [03]. Although the conventional CMSSM scenario requires more severe tuning than the 
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Sparticles 

Mass [GeV] 

(S)Particles 

Mass [GeV] 

m &. 

1682 

m c. 

2834 

m &. 

1530 


1157 

m o. 

581 

m i. 

2390 

rn u, 

1157 

m £ , 

2298 


1698 


414.5 

m u , 

799 

m No 

414.5 

m v , 

1636 

m N , 

414.5 

m V , 

1698 

M Ru 

1100 

% 

2298 

M n d 

298.5 


1682 

m 3 

550 

m Co 

1396 




Table 3: The soft scalar masses m/s , the up- and down-type Higgs masses M Rv d , and the 
gluino mass M 3 at the GUT scale. The subscripts Q a denote the mass eigenvalues for the 
left-handed Qi, up-type right-handed Ui, down-type right-handed T> t squarks, left-handed 
right-handed £ t charged sleptons, and right-handed sneutrinos N t with the three-generation 
i = 1,2,3. 


degree of 0.1 %, as pointed out in Ref. [2], certain ratios of the nonuniversal gaugino masses at 
the GUT scale relax the degree of tuning and observed 125 GeV Higgs boson mass at the same 
time. 


4.2 Relic abundance of the gravitino 

We are now ready to estimate the relic abundance of the gravitino. As stated in Sec. 13.51 
there are no significant gravitino productions from the inflaton, moduli, stabilizer, and SUSY- 
breaking fields after the inflation. However, there are two processes to produce the gravitinos 
associated with the decay of other particles in the MSSM. 

One of them is the decay from the thermal bath which is constituted of the relativistic 
particles after the reheating process. On the thermal bath, the dominant decay process comes 
from gauginos into gravitinos, because the couplings between the gravitino and other sparticles 
are more suppressed than those of gauginos as discussed in Refs. [331 05] • The abundance of 
the gravitino is estimated as 


nTP 

“3/2 


fc2 = E i + 


CL— 1 


M a {T R f 

3ml /2 


wMT r ? In (io^ev) ' 




(79) 


where w a and k a are the parameters whose values are defined in Ref. [45] and h is a dimensionless 
Hubble parameter. The thermal production of the gravitino is drawn in Fig. [2] in terms of the 
ratios of gaugino masses at the GUT scale Mqut, t'\ = Mi(Mqut)/M 3 (Mgut) and r 2 = 
M 2 (Mqu T )/M 3 (Mqu T ) with M 3 (M G ut) = 550 GeV. The gaugino masses at the reheating 
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r i 

Figure 2: Contours of the thermal abundance of the gravitino in the (r'i, r 2 )-plane. 

temperature M a (T R ) can be expressed as M a (M GU t) by employing the one-loop RG equations 
in the MSSM. The Planck Collaboration reported that the abundance of dark matter resides in 
the range of 0.1175 < O^/i 2 < 0.1219 [22] . where the upper and lower limits correspond to the 
dotted curves in Fig. [2l Here we assume that the dark matter only consists of the abundance 
of the thermally produced gravitino. 


NNLSP (Higgsino-like neutralino) 

Mass[GeV] 

Xs 

441 

NLSPs(right-handed sneutrinos) 

mass[GeV] 

R..2 

415 

2 

415 

U T 2 

415 

LSP (gravitino) 

mass[GeV] 

'h 3/2 

395 


Table 4: The masses of NNLSP, NLSPs, and the gravitino at the EW scale for the reference 
point (rq,^) = (6,3.5). The subscripts denote the mass eigenvalues for the sneutrinos (h), the 
Higgsino-like neutralino (y). 


The other process is the nonthermal gravitino productions from the NLSP and/or next- 
to-next-to-lightest supersymmetric particle (NNLSP). As shown in Table 0] when we take the 
ratios of gaugino masses (ri,r 2 ) = (6,3.5) consistent with the observed relic abundance of 
dark matter in Fig. [2] the NLSPs and NNLSP correspond to the degenerated sneutrinos and 
Higgsino-like neutralino, respectively. The relevant sparticle spectra are obtained by employing 
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the one-loop RG equations in the MSSM from GUT to EW scale with (ri,r 2 ) = (6,3.5) and 
the input parameters in Table [3] The full sparticle spectra are shown in the next section. Note 
that the degenerated sneutrinos do not have sizable interactions with the other (s)particles due 
to the tiny Yukawa couplings of Dirac-type neutrinos and then the soft SUSY-breaking masses 
of right-handed sneutrinos do not receive significant loop corrections. 

Since the gravitino and right-handed sneutrinos are weakly coupled with the other (s)particles, 
they are not thermalized. Thus, the nonthermal gravitino productions from the higgsino-like 
neutralino and sneutrinos are roughly estimated as 


a MTP h 2 = h 2 


(80) 


where m^o and are the mass and the thermal abundance of the Higgsino-like neutralino 
y 3 , respectively. The thermal abundance of the Higgsino-like neutralino is known to be small 
when the /i-term is smaller than wino and bino masses. Since the chargino and Higgsino- 
like neutralino are degenerated, both decay into the particles of the SM at almost the same 
decoupled time, which leads to the smallness of the thermal abundance of x 3 . After all, the 
nonthermal abundance of the gravitino can be neglected, 


n 


NT Pi 2 
3/2 n 


< 0 . 11 , 


(81) 


and the total relic abundance of the gravitino is approximated by the thermal abundance of 

itfl 


D 3 / 2 /r 2 — (82) 

However, the decays of neutralino and sneutrinos into the gravitino dark matter would 
threaten to spoil the successful BBN. The produced right-handed neutrinos via the sneutrino 
decay into the gravitino are suppressed due to the thermal abundance of y 3 , and then they 
are harmless for the BBN. On the other hand, the Higgsino-like neutralino decay into the 
gravitino affects the BBN. The authors of Ref. 03 suggest a way to relax the constraints from 
the BBN by assuming that the NLSP is the Dirac-type right-handed sneutrino. Although they 
consider the bino-like neutralino as the NNLSP, the sparticle spectra are almost the same as 
our obtained one. Because of the small thermal abundance of the Higgsino-like neutralino, it 
is then expected that our spectra are consistent with the BBN. 

Note that the nonthermal production of the gravitino is enhanced when the bino-like neu¬ 
tralino is NNLSP which corresponds to the small value of jr*! | in Fig. [2j In this case, it 
would break the successful BBN because of the large thermal abundance of the bino-like neu¬ 
tralino mm m- 


4.3 The Higgs boson mass, gravitino dark matter and sparticle spec¬ 
tra 

The ratios of gaugino masses at the GUT scale, rq and r 2 , are severely constrained by the relic 
abundance of the gravitino as can be seen in Fig. [2] In this section, we show that the mass 

8 In this paper, we do not take the gravitino production by the primordial black hole into account [46| . 
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Figure 3: The Higgs boson mass, the degree of tuning a /i- term, |A^| x 100(%), and the relic 
abundance of the gravitino 0 3 / 2 h 2 on the (rq, r 2 )-plane. In the blue shaded regions, the Higgs 
boson mass resides in the allowed range, 124.4 < my,, < 126.8 GeV [lj. The green dashed and 
solid lines show the 1% and 10% tuning, respectively. The red dashed curves show the relic 
abundance of the gravitino within ranges 0.1179 < Q 3 / 2 h 2 < 0.1215, reported by the Planck 
Collaboration [22] , 


of the Higgs boson further constrains the ratios of gaugino masses, rq and r 2 . The lightest 
CP-even Higgs boson corresponds to the SM-likc Higgs in the framework of MSSM. Without 
the loop corrections, the Higgs boson mass is much lower than the observed mass of the Higgs 
reported by the LHC experiment |48j. Although, the high-scale SUSY-breaking scenario is a 
simple solution as one of the possibilities to raise the Higgs mass, it requires the tuning to 
obtain the EW vacuum. Therefore, we consider the maximal mixing of left- and right-handed 
top squarks to raise the Higgs boson mass without a severe fine-tuning. 

With an approximation that the mass eigenstates of top squarks are nearly degenerate, the 
mass of the lightest CP-even Higgs boson cannot be realized. Thus, as pointed out in [2], we 
add the contribution from the mass differences between left- and right-handed top squarks in 
order to realize the observed Higgs boson mass and relax the degree of tuning at the same time. 
By employing the full one-loop RG equations of the MSSM from the GUT to the EW scale, we 
numerically calculate the Higgs boson mass which resides in the range of 124.4 < rrih < 126.8 [1 j, 
which is represented as the blue colored region in Fig. [3] and the degree of tuning a //-term, 
|A m | x 100%, is also given by the green dashed (1%) and solid lines (10%), respectively. From 
Fig. [3] there are the parameter spaces which are consistent with the relic abundance of the 
gravitino and the Higgs boson mass reported by the current cosmological observations [22j as 
well as the collider experiments |T] without a severe fine-tuning. 

In particular, at the reference point (rq, r 2 ) = (6, 3.5), the sparticle spectra, the Higgs boson 
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mass rrih, and the degree of tuning a /i- term |A M | x 100(%) are summarized in Tables 0 [5] 
and El It is then satished by all the experimental lower bounds from the LHC experiments 
for the masses of all sparticles in Refs. [13] and [49] . In general, the SUSY flavor violations 
are dangerous in the gravity mediated SUSY-breaking scenario due to the flavor dependent 
interactions. In our setup, there are vanishing H-terms and no flavor dependent soft SUSY- 
breaking terms at the GUT scale because the moduli do not have the F-terms. Even if the 
moduli have the F-terms at the vacuum, the SUSY flavor violations can be suppressed from 
the structure of the U(l)p charge assignments [50]• Thus, there are no serious SUSY flavor 
violations; especially, the decay rates such as n —* ey and b —* sy evade the present limits m 
152]. 
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Mass[GeV] 

Sparticles 

Mass[GeV] 
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2618 

h 

3241 

U 2 

2359 

e 2 
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Cl 
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2421 

C2 

2011 

/f 2 

2331 
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1735 

Tl 

2133 

h 
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t 2 

1447 

cl\ 

2625 

v ei 

3240 

d 2 

2620 

Y.2 
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Sl 

2522 


2330 

s 2 

2189 

^2 

415 

bi 

2117 

f'n 

2132 

^2 

1724 

v T2 

415 

Xi 

1723 

Xi 

444 

X.2 

1135 

X2 

1723 

X* 

CO o 

448 



Xa 

441 




Table 5: A typical sparticle spectra at the EW scale for the reference point, (r 1 ,r 2 ) = (6,3.5). 
The subscripts denote the mass eigenvalues for the following: up (h), charm (c), top ( t ), down 
(d), strange (s), bottom (6) squarks, the scalar electron (e), muon (p,), tauon (r), neutrino (z>), 
the neutralino (y), and the chargino (y 1 * 1 ). 
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m h [ GeV] 

mn[ GeV] 

mi [GeV] 

mn±[GeV] 

125.4 

1423 

1423 

1425 

A" 1 x 100(%) 

M^mzHGeV] 

M 2 {m z ){ GeV] 

M 3 (m z )[GeV] 

2.1 

1133 

1719 

1575 


Table 6: The neutral and charged Higgs boson masses and m//,, the degree of 

tuning a //-term, |A M | x 100(%), and the gaugino masses at the EW scale for the reference point 
( r i> r 2 ) = (6, 3.5). 
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5 Conclusion 


In this paper, we proposed the gravitino dark matter in the gravity mediated SUSY-breaking 
scenario based on the 4D AT = 1 SUGRA. The nontrivial Kahler metric of the SUSY-breaking 
held induces the mass hierarchies between the gravitino and the other sparticles for any value 
of the F-terrn of the SUSY-breaking held. Especially, the small Kahler metric of the SUSY- 
breaking held leads to the stable gravitino of mass (9(100) GeV with TeV scale gauginos and 
sparticles which would be the typical features in the natural MSSM with gravity mediation, if 
the gauge kinetic functions and the kinetic terms of the matter holds satisfy certain conditions. 
(See Ref. [21] for the case of CMSSM.) In the stable gravitino scenario, one can consider the 
low-scale SUSY without the cosmological gravitino problem, only if the NLSP decays do not 
spoil the success of BBN. 

As a concrete model, we considered the 5D SUGRA model on S 1 /Z 2 . Since the successful 
inhation mechanism as well as the moduli stabilization have been realized in 5D SUGRA [T9ll22j. 
we have estimated the moduli and inhaton decays into the gravitino dark matter. Although 
the produced gravitinos via the moduli decays seem to be dangerous from the cosmological 
point of view, their decays can be suppressed only if the moduli do not have the F-terms. 
Such a situation can be applied in our model, because the moduli, inflaton, and stabilizer 
fields have supersymmetric masses at the vacuum. Even if the supersymmetry is broken in 
the SUSY-breaking sector, their F-terms are suppressed by the gravitino mass at the SUSY- 
breaking minimum. When the NLSP and NNLSP are taken as the sneutrino and Higgsino-like 
neutralino, the nonthermal productions of the gravitino are negligible due to the small thermal 
abundance of the Higgsino-like neutralino. The smallness of the thermal abundance of NNLSP 
also relaxes the constraints from the BBN mmm, and at the same time, the amount of 
neutrinos via the sneutrino decay can be suppressed. Thus, the total relic abundance of the 
gravitino is approximated by the thermal abundance of it which depends on the gaugino masses. 
As pointed out in [2], the certain ratios of gaugino masses are also important to raise the Higgs 
boson mass in the MSSM without a severe fine-tuning. From Fig. [3] it is found that certain 
ratios of gaugino masses are consistent with the relic abundance of the gravitino as well as the 
Higgs boson mass reported by the recent Planck and LHC data [1, [22] , 

In this paper, we focus on the 5D SUGRA in order to show the realistic gravitino dark 
matter in the gravity mediation, and then the suppressed Kahler metric of the SUSY-breaking 
held is important to generate the mass hierarchies between the gravitino and other sparticles. 
When the 5D SUGRA is derived as the effective theory of superstring theories on a warped 
throat and/or M-theory on the Calabi-Yau manifold [33], the SUSY-breaking sector would be 
constructed from the gauge theory living on Dp-branes and/or NS5-branes. Especially, in the 
type II string, the visible and hidden sectors can be realized on the different D-branes which 
wrap the certain cycles in the internal manifold. In such cases, the different volumes of the 
internal cycles lead to the hierarchical Kahler metric between the SUSY-breaking held and 
matter holds in the visible sector. 
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A The F-terms of fields at the vacuum 


In this appendix, we derive the F-terms of the moduli, stabilizer, and SUSY-breaking fields 
at the vacuum by employing the reference point method. As discussed in Sec. 13.31 when we 
expand the fields around the reference points given by Eqs. (fT6j) and (l2T]h (j) —> </>| ref + S(p, 
(j) = T 1 ,Hi, X with I',i = 1,2,3, the Kaliler metric with Kahler potentials (IT4|) and ffl8|) are 
expanded by 


K — _i_ zyd) 

K u ~ K U + K IJ ’ 


where 
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in the field basis (T 1 , T 2 , T 3 , Hi, H 2 , H 3 , X), with 
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and the inverse of the Kahler metric is given by 


where 


K rI ' fI ' ~ (2ReT 7 ') 2 + 8ReT r Re5T r , K r Hi ~ -A h .8H u 
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A#; = (2ReT 


i/'\2 a ^i 


z , 4r = (2ReT / ') 2 ^. 


( 88 ) 


Here and hereafter, we omit the subscript of 0 at the reference point, that is, 0 = 0| re f- From 
the relevant expansions in the scalar potential (JTJ) with the Kahler and superpotential 

(TT5lh and (1IH1) , 

Drpj/\R 
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we obtain the scalar potential at the second order S(j), 

V ~ - 2 wW x (SX + SX) - 2^(2ReT / )wW T v H .{8Hi + SHi) 
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Finally, the extremal conditions for the relevant helds lead to the following variations of them: 
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and their F-terms become 
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, (i = 1,2,3), D X W ~ v. 


We also numerically checked these results, and then their F-terms can be suppressed by the 
tiny mass of the gravitino. 
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B The minima of fields during the inflation 

By contrast, the minima of fields during the inflation are different from those at the true vacuum. 
In this section, we derive the minima of fields by employing the reference point method. The 
reference points of fields during the inflation are chosen in the same way as those at the vacuum. 
Similarly, we expand the fields except for the inflaton Re T 3 around the reference points given 
by Eqs. (TT6l) and (l2Tj) . 0 —» 0| re f + 50, 0 = T 1 ', ImT 3 , Hi, X with I' = 1,2, i = 1,2,3. It is 
then supposed that H 3 is fixed at the origin due to the Hubble-induced mass. From the scalar 
potential (JT]) with the Kahler and superpotentials (CE3J, (EES]), and (ETHjl given by the following 
expansions: 

DrpjfW UJ ~\~ Wrpp H 5Hi X. rp J / rp J / W (ST 1 ' + Sf 1 ') + K t p(W H3 SH 3 + W x SX ) 

3 

+ K TllfI ,W H3 (8T T ' + ST I ')SH 3 + W T v T v H .5T r SHi + K T v ^ W T j' Hj 5T J '6H j , 

J'=j 

D Hl W ~ Wtih.ST 1 + K HlSl wSH 1 + K HiSi Wh 3 8HiSH 3 , 

D H2 W ~ W T 2h 2 ST 2 + K H2S2 wSH 2 + K H2B2 W H3 5H 2 5H 3 , 

Dh 3 W ~ W H3 + W t3Hs ST 3 + K H3 s 3 WSH 3 + d T s(K H3 s 3 )w(5T 3 + 8T 3 )5H 3 , 

_ |/y t3t3 T-J 

+ K H3 s 3 (Wh 3 \8H 3 \ 2 + W x 5H 3 8X) + 2 3 (ST 3 ) 2 , 

D X W ~ W x + KxxwSX + Ux(K X xM2\SX\ 2 + (5X) 2 + (5X) 2 ) + K X xW H3 5H 3 8X, 


I'=i 


W ~ w + WxSX + Wh 3 SH 3 + W t p n.ST 1 'SHi, 

3 

k -J2 


v=i 


v R eST 1 ' 1 / ReST J \ \ , l9 , x 

lllReT " + 2 () I + Z «‘\ SH <\ + Z *\ SX \> ( 94 > 


we obtain the extremal conditions for the relevant fields, and then their variations become 
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where 


5T 1 = dcr 1 + iSr 1 , SH, = dhi + i8k i: 8X = 8x + iSy, (96) 


with i — 1, 2, 3. 
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